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Abstract 

We consider Higgs pair production at the Large Hadron Collider (LHC) in a Randall- 
Sundrum scenario containing a radion field. It is shown that the enhanced effective coupling 
of the radion to gluons, together with contributions from a low-lying radion pole, can provide 
larger event rates compared to most new physics possibilities considered so far. We present 
the results for both an intermediate mass Higgs and a heavy Higgs, with a detailed discussion 
of the background elimination procedure. 
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1 Introduction 



The standard electroweak theory still awaits the discovery of the Higgs boson. After the 
Large Electron Positron (LEP) collider has set a lower limit of about 114.5 GeV on its mass 
[1], the responsibility of finding the Higgs now rests mostly on the Lage Hadron Collider 
(LHC). At the same time, puzzles such as the naturalness problem make a strong case 
for physics beyond the standard model (SM), just around or above the mass scale where 
the Higgs boson is expected to be found. It is therefore of supreme interest to see if the 
collider signals of the Higgs contain some imprint of new physics. This necessitates detailed 
quantitaive exploration of a variety of phenomenona linked to the production and decays of 
the Higgs. 

In this paper, we have studied pair production of the Higgs boson at the LHC as a 
possible channel for uncovering new physics effects. In particular, we show that the rate 
of such production receives a large enhancement in a class of theories with extra compact 
dimensions, namely, the Randall-Sundrum (RS) type of models containing a radion field. 

As has been mentioned above, the large hierarchy between the electroweak scale M EW 
and the Planck scale M P \ is somewhat puzzling. Whereas theories like supersymmetry 
and technicolour, each with its own phenomenological implications and constraints, have 
addressed this question, theories with extra spatial dimensions proposed as a resolution of 
this problem have recently drawn a lot of attention [2]. Basically, such theories postulate 
that all hitherto known particles and their standard interactions are confined to our familiar 
(3 + 1) dimensional spacetime (on a 'brane'), while gravity propagates in the 'bulk' including 
additional spacelike compact dimensions. On compactification, the (3 + 1) dimensional 
projection of gravity gives rise to both massless and massive graviton modes, where the 
latter interacts with the standard particles with a strength that can be perceptible in TeV- 
scale accelerator experiments. 

There are two major variants of the above approach. The first one of them [3] proposes 
a factorizable metric, large compact extra dimensions, and a compatification scheme that 
gives rise to a continuum of gravitonic modes on the Brane. The integrated effects of this 
continuum in observable processes have been the subject of a large number of phenomeno- 
logical studies. However, given the fact that we aim to have a bulk gravity scale in the TeV 
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range as a natural cut-off to the SM (and therefore as a solution to the naturalness problem), 
such an approach fails to explain a fresh hierarchy that it introduces, namely, one between 
the higher dimensional gravity scale and the inverse of the large compactification radius. 

The other approach is the one proposed by Randall and Sundrum [4]. In this model 
the fifth dimension corresponds to an S 1 /Z 2 orbifold and the (4+1) dimensional world is 
described by the following 'warped' metric 

ds 2 = e- 2KR ^% u dx^dx l/ - R 2 c d6 2 (1) 

where K is the bulk curvature constant and R c is the radius of the extra dimension. The 
theory postulates two D 3 branes, one located at 9 = where gravity peaks and the other at 
9 = it where the SM fields reside. The factor e ~ 2KR ^\ appearing in the metric is known as 
the warp factor. 

The interesting feature of such a theory is that projections of the graviton on the brane 
at (9 = 7r) gives, apart from the usual zero mode, a discrete massive tower with spacing 
on the TeV scale. In addition, the coupling of these massive modes to the SM fields is 
suppressed not by the Planck mass M P but by a mass M P e~ KRc7T . In fact, it can be shown 
that all quantities with the dimension of mass undergo such an exponential suppression on 
this brane, thereby offering a tangible solution to the naturalness problem with KR C ~ 12, 
i.e. with a combination of masses that are not too widely separated. 

The length R c in this scenario is called the brane separation, and can be related to 
the vacuum expectation value (vev) of some modulus field T(x) which corresponds to the 
fluctuations of the metric over the background geometry given by R c . Replacing R c by T(x), 
we can rewrite the RS metric at the orbifold point 9 = n as 

ds 2 = g v ; v s dx»dx u - T{xfd9 2 (2) 

where g™J = e~ 27rftrT( ^7^ = (^j^) Vnv Here f 2 = 24 ^ 5 an d M 5 , the 5-dimensional Planck 
scale. 

One is thus left with a scalar field 0(x). This field is called the the radion field [5]. Of 
course, the modulus field has no potential to start with. Thus one needs to generate a stable 
vacuum for T(x) at R c , which in turn can give (f)(x) a non-zero vev. This is done in terms 
of the Goldberger-Wise mechanism [6], using a bulk scalar field with non-zero vev's at the 
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two branes, whereupon a potential for the modulus field is generated, and one ends up with 
a radion field whose mass (m^) and vev ((0)) are both within the TeV scale. In particular, 
the radion can be lighter than the other low-lying gravitonic degrees of freedom. Thus it can 
very well act as the first messenger of a scenario with compact extra dimensions, and reveal 
itself in collider experiments. Several studies on the observable implications of the radion 
are available in the literature [7], [8]. 

Here we wish to focus on the pair-production of Higgs bosons mediated by the radion 
at hadron colliders. The two features that can be instrumental in enhancing the signal 
in this channel are (a) the accessibility of the radion resonance for > 2m h, and (b) 
the relatively enhanced radion coupling with a pair of gluons at LHC energies. Before we 
come to the details of the predicted signal, however, let us briefly recapitulate the various 
interactions of the radion in a theory of the above kind. 

We list the relevant interactions of the radion with the SM fields in the next section. The 
general features of Higgs pair-production via radion are discussed in section 3. Sections 4 and 
5 contain discussions on the predicted signals for m^ < 2mw and m^ > 2mw respectively. 
We conclude in section 6. 



2 Effective interaction of radion with the SM fields 

Radion interactions with the SM fields on the TeV brane (i.e. the one located at 9 = it) are 
governed by 4-dimensional general covariance. The radion essentially couples to the trace of 
the energy-momentum tensor of the SM fields in the following manner: 

Cim = -QjT£(SM) (3) 

where (0) is the radion vev. There are phenomenological limits on the m < / > -(4>) parameter 
space, from which the lower bound on (0) can range from about the electroweak symmetry 
breaking scale to about a TeV, while can in principle be even lighter than mw- 
The trace of the energy-momentum tensor of the SM fields is given by 



2m 2 w W+W-» - m\ZpZ» 



+ (2m 2 h h 2 - d^h) + ... (4) 



The photon and the gluons couple to the radion via the usual top-loop diagrams; an added 
source of enhancement of the coupling is the trace anomaly term [9]. This term arises 
from the fact the radion couples to the trace of the energy-momentum tensor, which can be 
equated with the four-divergence of the dilatation current associated with scale invariance 
of the theory. Scale invariance is preserved at the tree level in sectors of the theory which 
are massless and have no dimensionful couplings. However, quantum corrections can break 
such invariance, thereby giving a value to the four-divergence, which is proportional to 
the relevant beta-function. This term augments the coupling of the radion to the trace of 
he energy-momentum tensor. In general, the contribution of this term to the interaction 
Lagrangian can be expressed as 

R (n \ 

(5) 



T^(SM) C 



2g a 



For gluons, f3 s (g s )/2g s = — [a s /87i] bQCD where &qcl> = 11 — 2n//3, rif being the number of 
quark flavours. On the whole, the effective <figg interaction is given by 

iCts^ab 



[bQCD + IQCD} (P1-P2V^ - P\vP2y) 



(6) 



where p\ and p 2 are the 4-momenta of the gluons G° and G\. The function Iqcd is dominated 
by the top quark loop, and is given by 

2 
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Iqcd 

where x t = ^f, and 

X(x t ) = 



2 + 3^/x7^T X(x t ) - 2(x t - l)X 2 (x t ) 



(7) 



-sin 



Ti + i In 



, x t > 1 
+ 



x t < 1 



(8) 



l-y/T=jT ti 

It is important to note that for x t > 1, rif — 5 and hence 6qcd = 23/3, while for x t < 1, 
nf = 6 and the corresponding 6qcd = 21/3. 



3 Higgs pair production at LHC 

The pair production of a neutral Higgs boson from gluon-gluon fusion at the LHC has already 
been studied in the context of the SM [10] as well as several of its extensions such as the 
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minimal supersymmetric standard model (MSSM) [11]. As has been already mentioned, 
there is an expectation of some trace of physics beyond the SM being found in signatures 
of the Higgs, and this can only be possible by studying Higgs production and decays in 
as many channels as possible. From this standpoint, the pair production of Higgs deserves 
attention, because the predicted rate of such production is very small in the SM [10], and 
thus any excess can be interpreted as the signature of new physics. The mediation of the 
heavier neutral Higgs in SUSY has been shown to be the source of some enhancement in a 
region of the parameter space. Theories with extra dimensions, too, have been studied in this 
context, both in the ADD and RS models [12], where it has been reported that the mediation 
of gravitons can boost the Higgs pair production rates. What we wish to emphasize, however, 
is the fact that the presence of a radion in the RS context is of particular significance here. 
This is because (a) whereas the graviton resonance in the RS scenario is usually at too high a 
mass to be significant, a relatively less massive radion can be within kinematic reach, and also 
(b) the enhancement of radion coupling to a pair of gluons via the trace anomaly term jacks 
up the contributions. We have computed the predicted rates for both the cases of m# < 2mw 
and ran > 2mw, and analyzed the viability of the resulting signals with appropriate event 
selection strategies. We have used the CTEQ4L parton distribution function, setting the 
renormalisation scale at the radion mass. We have checked that the predicted results are 
more or less unchanged if this scale is set, say, at the partonic subprocess centre-of-mass 
energy. 

Using Breit-Wigner approximation for the resonant production of a radion, one finds the 
following expression for the cross-section for pp — > <fi — > hh (where the dominant subprocess 
is 99 -> -> hh): 

cr h h(pp -> <t> -> hh) = J dxi J dx 2 £l (xi) 92(x 2 )cr hh (gg -> -> hh) (9) 

where, 



gg] s[s 2 + Am 2 h s + Ami] 
128 m*W P - ml)* + m'l*] 



a hh ( g3 ^^hh)= s ^: do) 



where s is centre-of-mass energy for the partonic subprocess. s = Xix 2 s, s(— 14 TeV) being 
the proton-proton centre-of-mass energy, and Xi,x 2 , the momentum fractions of the first and 
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second partons (gluons) respectively, and denote the mass and total decay width of 
the radion [13]. 

Away from the radion resonance, the term proportional to T| is of little consequence, 
and the total rate falls as l/(0) 4 if the radion vev is increased. This is because the decay 
width of the radion in any channel is proportional to l/(0) 2 . Near resonance, on the other 
hand, the contribution is dictated by the term with in the radion propagator. Given the 
above dependence of the decay width on (0), the net contribution near resonance becomes 
practically independent of the radion vev. 

4 Event selection strategy and results: rrih < 2mw 

The signal of Higgs pair production depends on the final states produced by Higgs decays. 
We shall broadly consider two mass ranges here, one corresponding to rrih < 2mw, and the 
other, to rrih > 2rriw- 

In the first case, the dominant decay channel for each Higgs is H — > 66, so that the 
signal consists of four b-jets. Normally, such events are beset with a huge background arising 
from the following sources: 

• The QCD production of (66) (66), 

• The QCD production of (bb)(jj) with the two non-6 jets misidentified as 6, 

• The production of two(66) pairs from different partonic collisions, arising either from 
the same pair of protons or from two different pairs. 

• The electroweak production of Z(bb) and W(bb). 

A very relevant discussion of these backgrounds and the cuts that can be used to eliminate 
them can be found in [14]. These cuts have been adapted to our context in this study. 
Basically, the contributions from mistagged b-jets and different partonic collisions from the 
same proton pair are relatively small, the 'pile-up' effects arising from different proton pairs 
can be also handled with the help of effective b-tagging. For this, the following b-tagging 
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efficiencies have been assumed in different transverse momentum (p T ) ranges[15]: 



e b = 0.6, for p T > 100 GeV 

= 0.1 + p T /(200GeV), for 40 GeV < p T < 100 GeV 

= 1.5p T /(100Gel/) - 0.3, for 25 GeV < p T < 40 GeV (11) 

where we also assume that 6-jets are tagged only in the pseudorapidity region r\ b < 2. 

With such tagging efficiencies, it is found that the irreducible backgrounds become man- 
ageable with a further set of event selection criteria, particularly if the Higgs mass is not in 
close vicinity of the Z-mass (a fact that is confirmed by LEP data). The exact cuts used in 
our analysis are listed below. 

• Reconstruction of the Higgs boson mass: Since we do not differentiate between 
6-and 6-jets, there will be four b in the final state and hence one has three possible ways 
to pair them up. We choose those pairs which correspond to the smallest invariant 
mass difference between them. The reconstructed Higgs boson mass is then defined by 
Mh = [M blyb2 + M b3jb4 ]/2. Keeping in mind the finite resolution of the reconstruction 
procedure, this invariant mass has been smeared by a Gaussian distribution of width 
a ss \/M h . The window for the reconstructed Higgs boson mass in which our search is 
being carried out is 

0.9m Km - 1.5a < M h < 0.9m hM + 1.5a (12) 

• The mass difference between the invariant masses of the two pairs are required to be 
below the width of Gaussian smear: 

5M h = \M hlM -M hM \ < 2a (13) 

• Large azimuthal angles between the two jets of each pair in the transverse plane are 
demanded: 

54> blb2 ,5(f) b3b4 > 1 (14) 
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m, (GeV) m, (GeV) 

h n 

(a) (b) 



Figure 1: The Signal cross-section(pb) for the process pp(gg) — > <p ~* hh — > 6666 against the 
Higgs mass for m+= 250, 375 and 500 GeV and ((f)) = 500, 1000 GeV. 

• Events where the difference between the angles Sfo^, 506 3 b 4 is small are retained: 

\5<p hlh2 - 5<f) b:ib4 \ < 1 (15) 

• The b-jet are subjected to the following p^-cuts: 

p T , min > M h /4; p T:max > M h /A + 2a (16) 

• In order to further utilise the hardness of the b-jets in eliminating backgrounds, a 
minimum value is imposed on the 46 invariant mass My,: 

M Ab > l.9M h - 3(7 (17) 

After applying all the above cuts, the cross-section for the process pp(gg) — > 4> — > hh — > 6666 
is obtained for different radion vevs (0), radion masses and Higgs masses rrih- In Figure 
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l(a,b) we have plotted this cross-section against m h (GeV) (~ m h>in ) for (0) = 500, 1000 
GeV and for = 250, 375 and 500 GeV. 

The first conclusion to draw from the figures is that there is a substantial enhancement of 
the total rates over what is predicted in the standard model as well as in the case of MSSM. 
In addition, it also exceeds the predicted rates in the ADD-and RS-type models when such 
models assume graviton mediation to be the only new effect. This is particularly evident 
from the fact that the numbers presented here are after all the cuts have been applied and 
the b-tagging efficiency has been folded in, which effectively causes well below one per cent 
of the signal events to survive. 

As has been mentioned before, such an enhancement has two main sources, namely (a) 
the availability of the radion resonance, and (b) the enhanced coupling at the gluon-gluon- 
radion vertex. It can also be noticed by comparing figures 1(a) and 1(b) that the rates for 
Higgs masses corresponding to the onset of resonance are almost independent of the radion 
vev. Since we are showing here only that range of Higgs masses where the decay into the 
bb is dominant, such peaking effect is prominent only for = 250 GeV; for the higher 
values of m^, becomes large enough for the branching fraction in the bb channel to drop 
drastically, before the peak can be reached. 

Significance contours: 

With an integrated luminosity of 100 fb -1 , the above rates indicate a rather impressive 
prospect of detecting pair-produced Higgs bosons if radion mediation is operative. To gauge 
the actual situation, however, one must also remember that the backgrounds are never 
totally eliminated, and thus one must examine how the signal fares compared to the surviving 
backgrounds. This is depicted in figures 2(a) and (b) with contour plots in where, S and 
B are the number of events corresponding to the signal and backgrounds respectively with 
the above luminosity. In calculating the backgrounds, we have taken both statistical and 
systematic effects into account, assuming the systematic uncertainty to be 2% of the total 
background and adding it in quadrature to the computed background itself. The contour 
plots corroborate our expectation that the signal really stands out over a large region of 
the parameter space. Thus the so-called intermediate mass range acquires a high degree of 
visibility if a low-lying radion is present in the theory. 
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Figure 2: Contour plots in rrih - plane corresponding to (0) = 500, 1000 GeV and ^= 
= 5, 10. 

5 Event rates for > 2mw 

Let us now consider the case of a somewhat heavier Higgs, which can decay into a pair of 
W's or Z's. The signal in such a case consists in AW or AZ final states, where the SM 
backgrounds are negligibly small. In Figures 3 and 4, we have plotted the cross sections for 
the processes pp(gg) — > 4> — > hh — > W + W~W + W~ and pp(gg) — > — > /i/i — > ZZZZ (for 
cases where the VF and the Z decay into electrons and muons with total branching ratios 
of about 0.2 and 0.06 respectively) against the Higgs mass m h for different and ((f)). In 
addition, an average detection efficiency of 75% per lepton has been assumed. It can be seen 
from the figures that the AZ final states are very unlikely to be seen at the LHC. The AW 
final state, however, is quite substantial, especially in view of the absence of backgrounds. 
One may thus expect to see events ranging in number from a few tens to several thousands, 
depending on the radion mass, so long as the Higgs mass is within approximately 200 GeV 
For a radion mass on the order of 375 GeV, a slight peaking behaviour can be seen around 
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Figure 3: The Signal cross-section (fb) for the process pp(gg) — > — > hh — > W + W W + W 
where each of the produced W decay leptonically (W — > lui, I — e, fx) against the Higgs mass 
for m^= 375 and 500 GeV and (0) = 500, 1000 GeV. 



rfifi ~ 180 GeV. For a heavier radion, however, the decay width becomes so large that 
the peak is washed out. In this case, the r^m^-term in the Breit-Wigner propagator gives 
a substantial contribution even where the subprocess centre-of-mass energy is considerably 
away from m|. Such a contribution is responsible for the lack of dependence of the radion 
vev, following arguments given earlier. 



6 Summary and conclusion: 

We have looked at the pair production of the Higgs boson enhanced by radion mediation 
at the LHC For irih < 2mw, 4fc final states are investigated. It is found that in spite of 
substantial backgrounds to start with, a careful event selection strategy can lead to a 
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Figure 4: The Signal cross-section (fb) for the process pp(gg) — > <fi — > hh — > ZZZZ where 
each of the produced Z decay leptonically (Z — > l + l~,l — e, /i) against the Higgs mass for 
m^= 575 and 500 GeV and (0) = 500, 1000 GeV. 

discovery potential at the 5a or even 10a level. For > 2mw, on the other hand, one has 
to depend on the AW final states which are almost background-free. There, too, one should 
be able to see anything between about 100 and 800 events so long as the radion vev is within 
a TeV and the Higgs mass lies within about 200 GeV. Thus the presence of a radion can 
boost the very conspicuous phenomenon of Higgs pair production, over a large region of the 
parameter space, including the entire range of Higgs mass favoured by precision electroweak 
data. 
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